Decoherence and dephasing in Kondo tunneling through double quantum dots 
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We describe the mechanism of charge-spin transformation in a double quantum dot (DQD) with 
even occupation, where a time dependent gate voltage Vg{t) is applied to one of its two valleys, 
whereas the other one is coupled to the source and drain electrodes. The Kondo tunneling regime 
under strong Coulomb blockade may be realized when the spin spectrum of the DQD is formed by 
the ground state spin triplet and two singlet excitations. Charge fluctuations induced by Vg{t) result 
in transitions within the spin multiplet characterized by the SO{5) dynamical symmetry group. In 
a weakly non-adiabatic regime the decoherence, dephasing and relaxation processes affect Kondo 
tunneling. Each of these processes is caused by a special type of dynamical gauge fluctuations, so 
that one may discriminate between the decoherence in the ground state of a DQD and dephasing 
at finite temperatures. 

PACS numbers: 73.23.Hk, 72.15.Qm, 73. 21. La, 73.63.-b 



I. INTRODUCTION 

Manipulating spins in quantum dots in order to achieve 
new ways of processing information has been a subject 
of intense research in the last decadei One of the ba- 
sic building blocks for quantum computing is a double 
quantum dot (DQD)^ where the concepts of entangle- 
ment and phase accumulation may be effectively mod- 
elled and studiedi^ It turns out that the key problem, 
which hampers the progress in implementing the quan- 
tum information storage and processing is the impossi- 
bility of perfect isolation of a complex quantum system 
from the environment, which results in the loss of quan- 
tum coherence4i The phase coherence of tunnel trans- 
port through a quantum dot can be effectively controlled 
in an Aharonov-Bohm geometry- and/or in the Kondo 
regime^ To understand the generic features of the deco- 
herence phenomenon, various model situations are con- 
sidered, where dephasing and relaxation processes are 
controlled by specific mechanisms of interaction with the 
external environment. Among these processes the in- 
teraction with a phonon bathi and an electron liquid^ 
should be regarded as typical examples. 

The aim of this paper is to study a specific mechanism 
of decoherence and dephasing, which stems from the vio- 
lation of dynamical symmetry in a DQD. This dynamical 
symmetry is an intrinsic property of any quantum dot 
with even occupation, where the low-energy spectrum of 
spin excitations is formed by singlet (S) and triplet (T) 
statesii2*ii Since the interaction with a reservoir (elec- 
trons in metallic leads) violates the spin conservation 
law, S/T transitions accompany the cotunneling through 
DQD and thus contribute to the tunneling transparency 
and, in particular, to the Kondo-type zero-bias anomaly 
in tunnel conductance. 

The physical system in which we will study decoher- 
ence and dephasing effects is a DQD occupied by two elec- 



trons in a T-shape geometry (TDQD), where only one of 
its two wells is in tunnel contact with the leads (see Fig. 
^. As it is establishedfi2*ii the electron states in this 
kind of quantum dot may be tuned in such a way that 
the low-energy spectrum consists of two singlets and one 
triplet. The symmetry of the pertinent quantum states is 
described by a non-compact group, that is, SO{5). Deco- 
herence and dephasing will be studied through the appli- 
cation of a time-dependent gate voltage to the side dot in 
a TDQD, and its influence on the resonance Kondo tun- 
neling. Usually, decoherence in Kondo tunneling arises 
due to non-equilibrium spin flip processes induced by 
an external time-dependent potential Ji Here we propose 
an alternative mechanism, which involves processes with 
charge transfer between the two wells of the TDQD. As 
a result, the singlet excitons (E) are involved in dephas- 
ing and decoherence processes. It will be shown that 
these processes arise due to dynamical gauge fluctuations, 
which accompany singlct-triplct transitions within vari- 
ous multiplets of the SO (5) group. 

The structure of the paper is as follows: The under- 
lying time-dependent Hamiltonian is introduced in sec- 
tion II, flrst as a generalized Anderson impurity model, 
and then in terms of Hubbard operators. Section III 
is devoted to the derivation of an effective spin Hamil- 
tonian (following an appropriate Schrieffer- Wolff trans- 
formation) and presentation of time dependent scaling 
equations for the relevant exchange constants. Fluctua- 
tion corrections to the scaling equations are thoroughly 
discussed in section IV. This section is the central one, 
as it exposes the physical origin of dephasing, decoher- 
ence and relaxation in the present model system. Since 
it might appear too technical, it is following by conclud- 
ing remarks in section V, where the main results are ex- 
plained verbally. Related mathematical topics such as 
necessary ingredients of Group Theory, manipulations of 
Vertex Corrections and some aspects of analytical con- 



tinuation as used in the calculations are explained in the 
following three Appendices. 



II. MODEL AND TIME-DEPENDENT 
HAMILTONIAN 

We consider an asymmetric DQD studied in Ref. [13, 
where only the left dot is coupled to the leads (Fig. ^ 
left panel). The capacitive energies in the left and right 
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FIG. 1: Left panel: T-shaped double quantum dot (TDQD); 
right panel: Ground state and lowest excitations of doubly 
occupied TDQD renormalized by interaction with leads. 

dots are different, Qr ^ Qi so that the Coulomb block- 
ade completely suppresses doubly occupied states in the 
right dot. The gate voltages Vg'"^ are applied separately 
to the left and right dots. The new ingredient here is 
the introduction of a small "trembling" (time dependent) 
component in the gate voltage of the right dot, so that 



yr 



vr{o) + vgit). 



Let us first recall the spectrum of the TDQD in the 
absence of the time dependent component. The constant 
parts of the gate voltages are included in the energy levels 
of the double dot, eir = £i.r 
are tuned in such a way that 



+ V;'^''(0). These voltages 



€i + €r < 2ei + Q; < 2er + Qr- 

Hereafter, the Fermi energy in the leads is chosen as the 
zero energy level. As excited states of the TDQD with 
two electrons in the right dot are excluded from the low- 
energy part of the spectrum, the low energy levels of the 
TDQD arsia 



Et ^ ei + Sr - Mt, 
Es=ei+er-2f]V-Ms, 
EE^2ei + 2pV + ME, 



(1) 



where V is the potential which acts between the dots, see 
Fig.l. The notations T, S, E respectively refer to triplet, 
singlet and exciton states of the TDQD. In the exciton 
state, the two electrons reside on the left dot. The above 
results are obtained in first order of the small parameter 
(i = V/{ei — Er + Qi) <C 1. Anticipating further analysis, 
we include in the low energy spectrum the level shifts 
M\^T,s,E which result from renormalization of the dot 
levels due to the tunnel contact with the leads (the so 



called "Haldane renormalization"ii^ see below). It was 
shown in Ref. [13 that My > Ms, so that the singlet 
and triplet levels may cross due to this renormalization. 
Such level crossing occurs provided the charge transfer 
exciton energy Ee is not very high. Here we study just 
this regime, so that taking into account the above level 
renormalization, the ground state of the TDQD in a con- 
tact with the leads is a triplet state Et ( Fig. ^ right 
panel) . 

Having in mind the above energy level scheme |^ for 
the TDQD with constant gate voltages Vg'^{0), we now 
consider the influence of the trembling potential Vg(t) 
on the dot spectrum. As usual, one should discriminate 
between slow and fast components of a temporal per- 
turbation and treat the former processes in an adiabatic 
approximation. In order to separate adiabatic and non- 
adiabatic parts of the trembling potential it is useful to 
introduce the spectral density A{Q) of the trembling po- 
tential Vg'. 



Va{t)=' 



/ dnA{n)t 



iQt 



(2) 



In the present study we are mainly interested in the 
nearly adiabatic regime, where the spectral density is 
concentrated in the frequency interval < f2 < JIq with 
r^Q <C AsT = Es — Et- In this case the contribution of 
the trembling potential may be treated by means of time- 
dependent scaling approach, while non-adiabatic correc- 
tions should be considered perturbatively. 

The Hamiltonian describing the system schematically 
shown in Fig. ^has the following form: 



H 



H. 



dot 



H. 



(3) 



Here the first term is the band Hamiltonian, which de- 
scribes electrons in the leads [source (s) and drain (d), 
respectively] , 



H, 



hand 



L^ 



£kbCi^i,^Ckb. 



(4) 



k, a are the wave vector and spin projection, respectively. 
The tunneling Hamiltonian involves only electrons in 
the left dot. 



Hf. 



wj2i' 



ia^la ■ 



H.C.), 



(5) 



kcr 



where the operators di^ correspond to electrons in the 
TDQD. In the tunneling part[S|we have already excluded 
one of the two channels from the tunneling Hamiltonian 
by introducing even and odd combinations of electron 
wave function in two leads. Then, only the even standing 
wave Cka = 7f (cfcso- + Ckda) enters Htunifi) (we consider 
a completely symmetric device with lead-dot tunneling 
constants Wg = Wd = W). 

It is useful at this point to carry out some manip- 
ulations which will turn the discussion more transpar- 
ent. In the first step, a canonical transformation is 



pcrformcdrii^ii^ which ehminates the trcmbhng poten- 
tial from the diagonal part of the dot Hamiltonian 



Hdot — Hio + HlrO + Hr 



(6) 



with 



Hm = £ini + Qin^ 



1 1 



HirO=VY,{d\jr.+H.C.), 



Hr = ErTlr + Qr^-r + Vg{Q)nr + Vg{t)nr = iJrO + Hr{t). 

Here rij = dj^dj„ is the density operator for the electrons 
in the left (j = I)) and right (j = r) dot. The required 
transformation is, 






(7) 



with Ui = cxp[— i$i(t)rir] where the phase $i(t) is given 

by 

^^i^) = \j dt'vg{t'). 

This phase may be rewritten in terms of the spectral 
density introduced in Eq. 12 

<Pi{t) ^ J dn^iin)e'''' = vg J dn^-e'""'. (8) 

Then, expanding the exponent in terms of the weak trem- 
bling potential, we come to the expression, 



Hdot = H^Ji\ + iV^i{t)Y,{dl,di„-H.c.) 



Hfl + Hdot{t) 



(9) 



In the next step, it proves to be more convenient to 
work in a representation |A), which diagonalizes the dot 
Hamiltonian, i.e. to rewrite it in terms of Hubbard oper- 
ators X^'^ = |A)(A'|. After the standard diagonalization 
procedure for the time independent part of the Hamilto- 
nian we havaiS 

^doi = E^AX^^ A=T^,^,s. (10) 



Recall again that T, S*, E stand for triplet, singlet and 
charge transfer exciton states of TDQD, respectively 
while /i = ±1,0 are the spin projections of the triplet 
state. The diagonal Hubbard operators are constrained 
by the condition 



^^^^ = 1. 



(11) 



In terms of Hubbard operators the tunneling Hamilto- 
nian (0) acquires the form 



h[1^W l^ K,X'v- + H.c, 



h,aa' 



(12) 



Here the index To-' stands for the electron in the right 
well with spin projection ct', which remains in the TDQD, 
when the electron l^ escapes from the left well to the lead. 
As was mentioned above, the lead state ka means the 
even combination of lead electrons with the wave vector 
k. All tunneling transitions are described by the same 
constant W . 

Writing Hdot{t) from Eq. ^ in the form Hdot{t) = 
iV^i{t)Sir where 

Sir=Y,{dljia-H.c.), 

a 

one finds 

Sir\S) = V2\E), Sir\E) = -V2|5), Sir\T) = 0. (13) 

It follows from these equations that 

Sir = V2{X^^ - X^^) = iAV2. (14) 

The new scalar operator A is one of the generators of the 
50(5) group, which characterizes the dynamical symme- 
try of a biased DQD (see Appendix A). Then, the com- 
mutation with the diagonal Hubbard operators yields 

[Sir,X^^] = ^A^/2, [5,„ X^^] = -iA^/2 (15) 

(Sir is an anti-Hermitian operator). As a result, the GA 
transformation— leads to 



Hdot{t) = -^V^i{t)A. 



(16) 



III. 



TIME-DEPENDENT POOR MAN'S 
SCALING 



Next, the Haldane- Anderson scaling approachi^ should 
be applied to Hdot + Hband + Hum in order to rescale 
the energies in the TDQD. Here Hdot ® contains the 
time-dependent component. Unlike the static case con- 
sidered in Ref. llOl we deal here with a TDQD whose 
Hamiltonian is non-diagonal in A, such that the off- 
diagonal elements are time-dependent. However, for the 
slow trembling potential with characteristic frequencies 
ria <S^ Aes = Ee — Es, one may treat this time- 
dependent term adiabatically at least in zero-order ap- 
proximation. This means that before turning to the 
Haldane RG procedure, one has to get rid of the non- 
diagonal SE mixing terms with the help of a second 
time-dependent canonical transformation. This transfor- 
mation is given by the matrix U2 = exp 0, which is found 
from the condition^ 



Hdot (t) 



[e,Hl"j] 



ih 



de 

'dt' 



(17) 



Within the adiabatic approximation, one may carry out 
this diagonahzation procedure at each moment t neglect- 
ing any retardation and omitting the term in the r.h.s. 
of this equation. Then the matrix is given by 



Q^^V2V^l{tl^^ES 



A 



ES 



+ X^^). 



(18) 



At this stage, we apply the "adiabatic" Haldane RG 
procedure, i.e. calculate the scaling trajectories of the 
energy states, which evolve with the reduction of the en- 
ergy scale in the metallic reservoir from the initial value 
Dq to the actual value D. Since the transformations Q 
and (|17|) do not involve the triplet state, the scaling tra- 
jectory Et as a function of scaling variable ^ = Iu^Dq/D) 
is the same as in our previous calculations ii£ As a result, 
the triplet state includes only the time-independent Hal- 
dane shift, whereas the two singlet states acquire time 
dependence. 



Et ^ Er-MriO, 
Es{t) ^ Es-MsiO-Cs^lit), 
EE{t) -^ EE + ME{i) + Cs^l{t). 



(19) 



Here Cs = 2V'^ / /S.es is an additional time-dependent 
adiabatic part within the Haldane renormalization 
scheme. Thus, the term H^J^ in the dot Hamiltonian 



© is given by Eq. ^^ and 



Hdot{t) = Csn{t){X^''-X 



^SS\ 



(20) 



There is no correction to the matrix O from the time 
derivative in the r.h.s. of Eq. (|17|l at least up to first or- 
der in the adiabatic parameter ki = Vg{t)/ /S.es- Besides, 
due to the same transformation t/2, the components of 
the tunnel Hamiltonian containing operators X^'^ and 
X^^ become time dependent as well. 



Ht: 



Htim + Htun{t)- 



(21) 



The "trembling tunneling" contribution in these terms 
has the form 



Htun (t) 



Aes t" 



Sr^ 



X^^-)ck.-H.c. 

(22) 

Thus, the trembling gate voltage generates time- 
dependent tunneling through the singlet states, but the 
operator form of Htun is conserved. 

We see that the trembling potential involves a contri- 
bution of the excited singlet in the tunneling Hamilto- 
nian, and the latter introduces fiuctuations in tunneling 
through the ground state singlet. The spin S = I states 
are not involved in these trembling processes. However, 
eventually we expect that the triplet is subject to dephas- 
ing via the operators P and M (see equation 123 below) 
due to dynamical SO{5) symmetry of the TDQD (see 
Appendix A). One should remember that the state \E) 
enters the manifold of eigenstates of the isolated TDQD 



(I10|l . In the static case this high energy charge excita- 
tion is admixed with the low- lying singlet spin state | S) . 
This admixture changes the lead-dot tunneling rate and 
can lead to the inequality Mr > Ms mentioned earlier. 
We consider the situation well beyond the S/T crossover 
so that the positive singlet-triplet energy gap satisfies the 
inequality Ast = Eg — Ex ^ ^a- The adiabatic Haldane 
procedure ()19|l results in a time-dependent S/T energy 
gap 



AsT{t)=Es-ET~5{t), 
Sit) = Cs$?(i). 



(23) 



The adiabatic contribution of a charge transfer exci- 
ton (E) to the effective indirect exchange (cotunneling) 
Hamiltonian may be taken into account by means of a 
time dependent Schrieffer- Wolff transformation, which 
includes Htun{t) l|22|) . This transformation is described 
in Ref. |3 for the dot occupied by a single electron with 
S — 1/2. In our case the exchange Hamiltonian is af- 
fected by the trembling perturbation in spite of the fact 
that the localized spin S = 1 is not directly affected by 
the time dependent potential. However, charge fluctua- 
tions induced by time dependent tunneling H22|l perturb 
the S- and ii^-states, which are connected with the spin 
by the kinematics of vector and scalar operators of the 
dynamical symmetry group 5*0(5). All these corrections 
can be incorporated in the Schrieffer- Wolff (SW) Hamil- 
tonian Hsw for the dot obeying the SO (5) symmetry. 
The time-independent part of the SW Hamiltonian takes 
the formic 

H 



(0) 
SW - 


tt(0) rjW 


(24) 


r(0) 


= i (^tS^ + EsV^ + EeM^) + Q{N - 


-2)^ 


)) 

}tU7l 


= J^S ■ s + Jo^^P • s + Jo^^M • s. 





Here the terms describing the irrelevant potential scat- 
tering are omitted (see, however, the discussion of fluc- 
tuation corrections in Section IV). 

The Hamiltonian H24|l is expressed in terms of SO (5) 
group generators presented in Appendix A. Here N is the 
operator of particle occupation number, and the last term 
in H24|) describes the constraint N = 2 imposed by the 
Coulomb blockade Q on the charge sector of the Hilbert 
space. As we mentioned already, the operator S is not af- 
fected by the time-dependent canonical transformations. 
However, the operators P and M will be involved in the 
time-dependent processes. 

We generalize the adiabatic SW transformation de- 
rived in Ref. |3 by introducing the matrix C/3 = exp iT 
with T defined as 

T(t) = E Hit)X'''ck. + v^{t)X^^^cu. - H.c] . 

k(7 

(25) 
The coefficients v^ are chosen to eliminate the time de- 
pendent terms H22|l in the effective Hamiltonian 

Hsw == cyi]i{-iT)[Hdot + Hband + Htun] cxp iT, (26) 



where the first and third terms in the Hamiltonian de- 
pend on time. To fulfil this request, the transformation 
matrix should satisfy the equation 



Htun + [Tf, {Hdot + Hband)] = ^^ 



dT 
'dt 



(27) 



(cf. Eq. I17|l . Repeating the procedure used in the calcu- 
lation of the matrix U2 and including the first order cor- 
rections in the adiabaticity parameters ref (i) = Vg{t)/ei 
and Avf(i) = Vg{t)/{2ei — e^) from the time derivative on 
the r.h.s. of Eq. (|27|) . one finds the following expression 
for T{t) 



T{t) = TT + Tsit) + TEit), 

W 
T- 



(28) 



T5(t) 



ei - Mt 

kaa' iy=0,l,l 

W - wl(t) 



ei-Ms{t) 



=0,1,1 

^(X^-'-^Cfc^-i/.c), 



ka 



The time-dependent quantities in T5 and T^ are 
Ms{t) = Ms - Cs^lit), Msit) =. Me + Cs$?(i) 



s ^^ _ VWV^'^iit) 



w^{t) 



^ES 



wiit) 



VWV2^f{t) 



^ES 



with $1 (i) = $1(0 - «f (i) and $f (i) = $1(0 - af (i). 
This procedure results in the appearance of a time- 
dependent component in Hcotun , which modifies the cou- 
pling constants J^^ a.ndJ^^ in the Hamiltonian (|24|l . As 
a result, we get the following effective exchange Hamilto- 
nian which induces the decoherence and dephasing effects 
in the Kondo tunneling through TDQD: 



Hr. 



Here 



Jo' ' S • s + J"" (i)P • s + J^' (t)M • s 



(29) 



j5T(^) _ W- wljt) Wjl - V/AEs)V2<^i{t) ^ ^3Q^ 



ei - Ms(t) 



^ES 



J^'it) = 



W-w^{t) W{l~V/l^Es)V2<^i{t) 



2er- ei+MEit) 



^ES 



(31) 

and the time-dependent parts of the coupling constant 
J^^ and J^'^ may be obtained by expansion of its general 
form in the adiabatic parameters ki, wf , k^. 

The first manifestation of this time dependence, which 
can be seen already within an adiabatic approximation 
is the uncertainty in the definition of the Kondo temper- 
ature. To describe this effect, we refer to the poor man's 



scaling equations for the Kondo effect in DQD for time- 
independent Hamiltonian l|24|) derived previously^ In 
the adiabatic regime, these equations retain their form, 
but the coupling constants depend on time t as a param- 
eter. In particular in the limit when the exciton state 
\E) is quenched, the system of scaling equations has the 
form 



dji/d( = -j!-Mt)]^; dj2/d( = -2jij2{t). (32) 

, Po is 



TT ., _ „„ J ST ,„ _ p^jET^ p^ ig ti^e 



Here ji = po^ , h = PoJ"^ , J3 
density of states at the Fermi level. 

As is already established in the theory of Kondo-effect 
at a singlet-triplet (S/T) crossingfiSiiSiiSiii in the limit 
AsT S> Tk the solution of these equations may be ex- 
pressed in terms of T^'o = -Do cxp[— l/(ji + J2)]: 



Tko 



T, 



A'O 



iC 



AsT(i) 



(33) 



where <^ < 1 is a constant. The relative amplitude STk 
of adiabatic variations of the "time-dependent Kondo 
temperature"— may be estimated from this equation: 



ST, 



K 



Tko 

Ast(O) 



i+C 



2Cs^l 



(34) 



(<t> J is the mean square value of the trembling parameter) . 
In this asymptotic region, the time-dependent corrections 
are insignificant. In accordance with the above men- 
tioned theory of an S/T transition for a time- independent 
gate voltagejiSii^iiSiii the Kondo temperature increases 
with decreasing Ast and reaches its maximum Tk ~ Tko 
at a crossing point Ast = 0. One may expect that the 
role of non-adiabatic corrections increases with decreas- 
ing Ast- As a result, the scaling behavior should be 
seriously violated when approaching the crossing point, 
and the deviations from the prescriptions of an adiabatic 
theory in the dependence Tk{Ast) should grow accord- 
ingly 



IV. FLUCTUATION CORRECTIONS TO RG 
EQUATIONS 

In this section we go beyond the adiabatic approxima- 
tion and take into account decoherence and dephasing 
corrections to the Kondo tunneling. Unlike the mech- 
anisms described in Ref. 0, where the time-dependent 
spin-flip processes were the source of dephasing, we ap- 
peal here to gauge fluctuations, which arise because two 
singlet states \S) and \E) are involved in the formation 
of a Kondo resonance in a triplet ground state \T) in the 
process of time-dependent cotunneling. 

The first source of non-adiabatic corrections is associ- 
ated with the fluctuations of the energy gap Ast iPSjl - 
which may be converted into gauge fluctuations of a 
Casimir operator ljA7p . Using an obvious chain of sub- 
sets, 

C/(l) C SU{2) c 5*0(4) c SO{5), 



we conclude that the trembUng potential affects the least 
continuous subgroup of 5*0(5). 

To perform this conversion we turn to the fermionic 
representation of the generators of the 5*0(5) group (see 
Appendix A). We now work in the subspace of low-energy 
states \T),\S), which are described by the dynamical 
symmetry group 5*0(4). Then the non-adiabatic fluctua- 
tions of the energy gap Aes raay be described as fluctua- 
tions of the "reduced" Casimir operator S^ 4- P^, which, 
in turn is rewritten as a local constraint (|A10() where 
the last term is eliminated. We introduce the "fluctu- 
ating constraint" in Hdot by means of a time- dependent 
Lagrange factor jjl: 



with 



H'^,,^Haom-^x{t){R-l) 



R = flfs+ J2 /-i^ 



(35) 



Thus, the fluctuating constraint takes into account time- 
dependent admixture of the i?-state to the S'-state. 
Although the exciton E is excluded from the effec- 
tive spin space, its effects arc present through a "ran- 
domization" of the constraint. The fluctuating time- 
dependent Lagrange factor can be eliminated from the 
5*0(4) fermionic Green's functions by means of the global 
U{1) gauge transformation performed in all sectors of the 
/-representation 



/a -^ fxe 



■■4{t) 



(36) 



Minimizing the free fermion part of the Lagrangian, one 
has h((){t) = fi{t) (here A = 1, 0, 1, s). 

On the other hand, the constraint fluctuations in terms 
of time-dependent Lagrange factors can be reformulated 
as a propagation of fermions along the time axis in the 
presence of a time-dependent external scalar potential. 
The solution of this problem results in the appearance 
of an imaginary part in the zero frequency spin fermion 
propagators, ImEj(O) = 7 (see below). Eventually, the 
existence of this damping prevents the formation of a 
coherent Kondo tunneling regime through DQD at zero 
temperature, so the effects connected with this type of 
gauge fluctuations can be qualified as decoherence. 

Another source of non-adiabatic gauge fluctuations 
comes from the time-dependent poor man's scaling so- 
lution. The co-tunneling Hamiltonian (|31() contains 
the time-dependent corrections to P • s and M • s. In 
this case (unlike the flrst mechanism), the source of 
gauge fluctuations is due to non-diagonal operators P, 
M, so that the time-dependent coupling constants are 
parametrized as J2(i) = .72 e**^^^*^ and J3(i) = jae**^^^*^. 

The phases in j^^ can be eliminated by the local U{1) 
gauge transformation 



f\ ^ f\e 



i-dx{t) 



(37) 



performed in the 5-, T- and £*-sectors of the fermionic 
representation of the 50(5) group. The local gauge 



phases may be represented as 

OsT{t)='l}T{i)--&s{t), eET{t) = Mt)~-SE{t) (38) 

SO that the phases i^a(0 have different time derivatives. 
Therefore, the effects related to this type of gauge fluc- 
tuations are, in fact, dephasing effects stemming out of 
weakly non-adiabatic TS and TE transitions. 

Now we turn to the calculation of gauge fluctuation 
corrections in the self energies and vertices entering the 
Kondo tunneling diagrams. To calculate the damping of 
retarded spin-fermion propagators in the triplet sector of 
phase space 

GT,{t t') = {f,{t)fl{t'))R = -i{[f,{t)fl{t')]+) (39) 

we apply perturbation theory for time-dependent exter- 
nal potentials to the last term in the Hamiltonian H35|l . 
In the corresponding diagrammatic technique the times t 
are marked by crosses and the correlation functions con- 
necting the fluctuations related to the instants i, t' are 
denoted by wavy lines. Then the first order diagram for 
the self energy of GTfj.{t~t') is given by the diagram Fig. 
121 (the spin-fermion propagator is represented by the full 
line). We consider the simplest case of white- noise cor- 
relations lT?{<f>i{t)(j)i{t')) = roS{t — t'). The correspond- 



j^-^i. 



-^ 



%f- 



t ^i t' 

FIG. 2: Self energy correction due to gauge fluctuations in 
spin-fermion propagator. 

ing contribution to the self energy is a purely imaginary 
damping 

71 = h/2Td ^ O|ro. (40) 

This first order approximation is valid as long as the 
damping is weak in comparison with the Kondo tem- 
perature, h/Td <?C Tk- If these two quantities are com- 
parable then the whole sequence of " rainbow" diagrams 
should be inserted in the spin-fermion self energy. This 
means that instead of the white noise we get a more re- 
alistic description of trembling, which takes the retarda- 
tion effects (non-gaussian corrections) into account. This 
specification is not too essential for our purposes. More 
important is the low-energy cutoff ~ h/rd in the RG pro- 
cedure, which prevents crossover to the strong-coupling 
limit of the Kondo tunneling due to incoherent phase 
fluctuations. 

Our next goal is to flnd non-adiabatic corrections for 
exchange vertices, which contribute to the dephasing 
mechanism mentioned above (sec Eq. I38II . Dephasing 
means that the gauge transformation (|37() eliminates the 
phase 9sT within the accuracy of phase fluctuations, i.e. 
the local gauge phases have the form 



9ETit)=Mt)-^Eit) + ^eit). 



(41) 



Then, expanding the exponents in the non-adiabatic ver- 
tex corrections J^^ [e^f^W - l] and J^^ [e'f^(t) - ij , 
we obtain the fluctuating part of the Hamiltonian H29(l 
in a form 



5Hc 



= i [j^'^ips{t)V ■ S + J^^(y9e(t)M • s] . (42) 



Adiabatic vertices H24() and non-adiabatic corrections 
(|42|l may be presented in a graphical form (Fig. ^. The 
bare interactions are shown in Figs. inia,b) by means of 
four-tail vertices where the operators S, P and M are 
expressed via fermion operators in accordance with Eqs. 
(|A9() . Here solid lines stand for spin fermions and dashed 
lines represent conduction electrons. 

The time-dependent vertices (|42|) are shown in Fig. 
O;, where the broken line stands for the fluctuation field 
D{t) = (^(t)^(0))fl. 





y \ 



a) 



y \ 



b) 




FIG. 3: (a) Adiabatic exchange vertices for triplet channel 
T; ( b) Similar diagrams for ST and ET channels; (c) non- 
adiabatic corrections to two latter vertices. 

If the fluctuating trembling signal is characterized by 
retardation, 

D{t) = -ia^g-Cltl (43) 

then the Fourier transform of this correlation function is 



THj 




(a) 



Tn, 



T^^l 




S 
(b) 



T^, 



FIG. 4: First leading parquet diagram (a) and first non- 
adiabatic correction to it (b). The pseudofermion, electron 
and fluctuation propagators are represented by solid, dotted 
and broken lines, respectively. The absence of direction in 
electron propagator lines assumes that both directions (clock- 
wise and anti-clockwise) are possible. 



tion of the diagram in Fig. ^p is parametrically small 
compared with that of Fig. 0^ due to a small factor 
{aj^'^Y ^ 1; this diagram represents a building block 
for construction of non-adiabatic corrections for vertices 
and self-energy parts. These corrections affect both adia- 
batic vertices (see Appendix B) and the self energy S]t((j-') 
of the spin- fermion propagator Gt^ (Fig. 0). The latter 
diagram may be obtained from the vertex (Fig^Jj) by 
gluing two free electron lines in the electron propagator. 
Besides, the diagram (Fig. [SJ is connected with first non- 
parquet TT vertex by a Ward identity (see Appendix B). 
We will use this identity below, when calculating the spin 
relaxation corrections. 




FIG. 5: First non-adiabatic correction to propagator Gth 



D{u) 






(44) 



Other elements entering the vertices are GTy{t -^ t') de- 
fined in Eq. (|39|l . a similar propagator for spinless 
fermions representing the singlet state 



G5(t-t') = (/.(i)/I(t')> 



B.- 



and conduction electron propagators 

gkAt-t')^{c.kait)cl^it'))R. 



(45) 



(46) 



The Fourier transforms of the Green functions (|39|l , H45|) 
and H46|l arc 

gk„{uj) ^ {u - Eka + iriV^ ■ (47) 

(here -q -^ +0). 

The first fluctuation correction to adiabatic vertices 
Fig. 2|a, is displayed in Fig. ^p. Although the contribu- 



The imaginary part of the self-energy Fig. [S] is given 
by the following expression 



tt{u)) 



ll\lYiT{tjj) = (ctJ2) X 



(48) 



deide2lmKi%iei - S2)liaD{uj - ei)lmGf\s2) 



I tanh -r:^ + coth — rrr— I I tanh -r:^ + coth 



2T 



2T 



2T 



2T 



We assume that the decrement 7 in the propagator D{t) 
(I43|l is small in comparison with the energies Ast and 
A^;^ since we remain in the weakly non-adiabatic regime. 
Here K^{llj) is the Fourier transform of the spin suscep- 
tibility of the electron gas in the leads 

X«(i-t')= J2 (4.,Wc....(i)cl3,Jt')cfc,.,(t'))fi 
which may be expressed via convolution of two electron 



propagators, 



/V «(t - t') = y" 9ka, {t ~ t')gk+q,a^ {f - t). 



kq 



The actual frequency interval is very narrow in compar- 
ison with the electron Fermi energy, oj ^ ep)- At low 
frequencies the Fourier transform of the electron-hole po- 
larization loop K^-{uj) behaves as 

Im K"{uj) ^ pq — . 

Sp 

Using this relation, one may estimate the dephasing rate: 



Asre"'^^^/^, T < A 



h 



ST 



72 = 



2tt{uj -^ 0) 



(aJ2 



T2 



iST 



AsT < T < A 



ET 



We see that this contribution to the dephasing effect is 
frozen out at low temperatures because the singlet states 
responsible for dephasing is depopulated at T ^ A^s- 

Now the self-energy corrections to the spin-fermion 
propagator Grfj. represented by the diagrams Fig. Inland 
Fig. Omay be used in the scaling equations (|32|l . As a 
result of fluctuation corrections, Grfi acquires the form 



GT^(a;) = (u + i^y 



(49) 



where 7 = 71 at T — > and 7 = 72 at T — > Ast- These 
propagators should be inserted in the diagrams arising 
in the poor man's scaling proceduroiS, (the first of these 
diagrams is shown in Fig. \^. Then wc immediately 
conclude that this imaginary part transforms into the 
infrared cutoff of Kondo singularity. This means that 
there is no significant effect of fluctuations on Kondo 
tunneling at high T ~ Ast provided 72 <C Tk-, where 
72 ^ T^/AsT- However both dephasing and decoher- 
ence prevent the achievement of the unitarity limit at 

Our last task is to clarify the contribution of vertex 
corrections presented in Figs. ^Jj, El and [7\ to the scal- 
ing equations H32|l . It is clear that the real parts of these 
corrections only slightly rcnormalizc the coupling param- 
eters ji^2 and do not influence the scaling trajectories. As 
to the imaginary part Tj,j,, it determines the longitudi- 
nal 1/Ti and transverse I/T2 relaxation rates by means 
of the correlation function, 



Xt(^) 



(S,S) 



(50) 



which can be interpreted as line-shape. While all dia- 
grams Fig. 6 a-i contribute both to l/Ti and l/r2 re- 
laxation rates, the diagram Fig. 7. a corresponds only to 
scattering processes without spin flip, resulting in slightly 
different temperature behavior of longitudinal and trans- 
verse relaxation rates. We note, however, that the dif- 
ference between 1/Ti and I/T2 appears only beyond the 
leading logarithmic approximation. 



The imaginary part Tgrp is associated with spin relax- 
ation processes determined by the specific kind of dy- 
namical spin susceptibility 



Xst{(^) 



(P,P) 



(51) 



which describes the response of a magnetic system with 
SO(n) symmetry not to external magnetic field but to 
the perturbations intermixing triplet and singlet compo- 
nents of the spin manifold [including gauge fluctuations 
()41|l ]. This correlator, given by the irreducible S/T loop 
(see Fig.7.b-d) leads to the appearance of a new I/T3 re- 
laxation rate associated with the inverse time of transi- 
tions between the singlet and the 3- fold degenerate triplet 
state. 

The diagrams (a,b,d,e,g,h) presented in the first two 
columns of Fig. are in fact taken into account in the 
RG equations together with the diagram (Fig. ^Jd). Each 
of the three diagrams (c,f,g) in the last column of Fig. 
contains two singlet propagators Gsic), so one should 
consider them together with the non-parquet diagram 
(Fig. [7^). This vertex correction to J^'^ together with 
non-logarithmic corrections (Fig. [Tja-d) due to fluctua- 
tion induced (T -^ S —^ T...) transitions accompany the 
adiabatic TT- and S'T-exchanges and introduce the lon- 
gitudinal and transversal spin relaxation times in the RG 
procedure (cf. Refs. ,,8 i|). 

The fluctuation induced vertex correction (Fig. [7^) 
may be estimated with the help of the Ward identity for 
Gs{(jJ -^ 0) written in the form 



aE 



TT 



dA 



TS 



rTT(o), 



(52) 



where T^tt is the self energy shown in Fig. and Ttt 
is the triplet vertex (Fig. [7|(a)). The Ward identity in 
this context corresponds to spin conservation in the pro- 
cess of quasi-elastic scattering. We conclude from these 
estimates that the non-adiabatic corrections ?i/(2ttt) = 
ImTxT presented by the diagram Fig. [71(a) are exponen- 
tially weak at low T <^ Ast similarly to the dephasing 
rate 72. 

Let us now consider the fluctuation corrections to the 
vertex J^^ represented in Fig. [7|b-d. There is no Ward 
identity for this inelastic process, which docs not conserve 
spin projection, so we calculate the vertex Tst straight- 
forwardly with the help of Eq. ljG5p . The imaginary part 
of the diagram Fig. [7}p gives the following correction to 
Tst{^,0) at w « Ast (see Appendix B): 

r T/Ast, Ast<T< Aet 
MPoTfr] = («'j1) <^ (53) 

[ uj/AsT, T < Ast 

The estimates of the next diagrams Fig. [7|;,d give a sim- 
ilar result. Such behavior of non-diagonal vertex correc- 
tions is predetermined by the threshold character of T/S 
excitations at finite frequency^ 



These vertex corrections should be inserted in the 
scahng equations H32|) . The imaginary part of the ex- 
change vertex introduces an additional cutoff in the scal- 
ing procedure;^ Taking into account the fact that the 
contribution of the singlet state to the flow equations 
controlled by the vertex j2 is frozen for D < Ast, one 
immediately finds that the relaxation processes practi- 
cally do not influence the cutoff of J2 because 

h/{2TsT) = ImTsTiuJ = Ast) < Ast- (54) 

The real part 

RcborsT] - («'j|) (T/Ast) HD/T), (55) 

just slightly disturbs the flow trajectories. Hence, it may 
also be neglected in perturbative estimates. 

As in the case of non-equilibrium Kondo effect^, the 
spin-relaxation processes are controlled by the imaginary 
part of the susceptibility K^'{lo) of the electron liquid 
in the leads. However, in our case the susceptibility is 
moderated by the gauge fluctuations (broken line inser- 
tions in the electron-hole or singlet-triplet loops in Fig. 
Efa-d)). As a result, the relaxation corrections do not af- 
fect Kondo tunneling at least in the weakly non-adiabatic 
regime. 



V. CONCLUSION 

A double quantum dot with a weak trembling poten- 
tial applied to the right dot (Fig. ^) turned out to be an 
excellent model system in which all facets of decoherence 
phenomenon are exposed as observable effects. From a 
theoretical point of view, this system is especially attrac- 
tive because decoherence, dephasing and relaxation are 
induced by the same gauge fluctuations, which develop 
in the constrained Hilbert space of the spin manifold 
{T, S, £'} (Fig. QJd) coupled to a Fermi bath of conduc- 
tion electrons. All these processes may be discussed in a 
general context of the theory of decoherence in quantum 
systems in contact with a thermal bath. 

The decoherence effect characterized by the time Td 
(I4(JII is related to the structure of the ground state of 
TDQD in contact with the Fermi bath. It may be in- 
terpreted in terms of the Superselection Rule introduced 
by Wick, Wightman and Wigner for the description of 
baryonic charge (see Ref. |23for a recent review). Indeed, 
there is no symmetry ban for superposition of two sin- 
glet states \S) and \E), but this superposition arises only 
as a result of the coupling of the TDQD with the Fermi 
bath. The covering group, which describes the symmetry 
of the manifold {T, S*, £'} is 5*0(5), and the dynamical su- 
perposition is controlled by [/(I) gauge fluctuations 13()|) 
under the Casimir constraint (|A10|) . Due to these fluc- 
tuations, the coherent Kondo- type ground state of the 
system 'TDQD + Fermi bath' cannot be reached. 

The dephasing effect characterized by the time tt 
stems from the phase averaging in thcrmodynamical en- 
semble at finite temperature. Dephasing 9sti Set l|38f) 



emerge in a process of exchange scattering induced by 
the random trembling potential. The scattering proba- 
bilities are added incoherently, so that the spin-fermion 
self energy acquires an imaginary part (Fig. \^. Such 
processes are generally classified as decoherence induced 
by dressing of bare statesiSi 

The relaxation effects characterized by the times ttt 
and TsT remind those known in the conventional theory 
of spin relaxation. Although we describe them in terms 
of triplet-triplet and triplet-singlet transitions, one may 
reinterpret the same processes in terms of longitudinal, 
transversal and S/T relaxation rates 1/Ti, I/T2 and I/T3 
because both TT and TS processes contain spin conserv- 
ing and spin reversal components [see Eqs. (|A9() ]. How- 
ever, since all these processes are controlled by the small 
coupling constant {aJ2 ) induced by the trembling poten- 
tial, the relaxation contribution to dephasing processes is 
ineffective. 

To conclude, we have shown in this paper that charge 
fiuctuations can be transformed into spin fluctuations, 
which result in decoherence and dephasing of the Kondo 
effect in the double quantum dot system. All these phe- 
nomena arise due to intrinsic dynamical 50(5) symme- 
try of the spin multiplet. Although we considered only 
a weakly non-adiabatic regime where decoherence effects 
are small by definition, the results are instructive, be- 
cause one may strictly discriminate between pure deco- 
herence of the ground state and the finite temperature 
dephasing and relaxation effects in a situation, where all 
these effects are due to gauge fiuctuations within a spin 
multiplet. Strongly non-adiabatic response to trembling 
gate potential demands special consideration. 
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APPENDIX A: LIE ALGEBRA FOR 
ASYMMETRIC DQD 

If the excitonic state is included in the set of the energy 
levels, then the transitions between different states are 
described by the 05 algebra. In addition to standard 
5 = 1 operators 



5+ = V2(X 



10 



X 



11 



X 



-1,-1 



5- = \/2 (A°i + A- 



(Al) 



10 



one should introduce two more vectors. The vector P and the Casimir operator 
with the spherical components 

P, = -(X°^ + X^"). (A2) 



- M^ + ^2 = 4. (A7) 



defines transitions between the singlet \S) and the com- It is important to remind once more that the Hamil- 

ponents \T^) of spin triplet. Similarly, the vector M with Ionian Hdot does not commute with vectors P and M 



components because 



[P„p2] = -[P„S2] = 2(X'5o - X°^), 



^^ = -(^ +^ )• (A3) [P+,p2]=-[P+,s2]=4(Xi^ + X«i) (A8) 

determines transitions between the left exciton and the 

triplet. The 05 algebra is completed by the operator A The fermionic representation of SO{5) groupiiS^ is 



A^-iiX''" ^X"""). (A4) 



characterized by 5-vector q^ = {fLifoJh flJl) 



The Lie algebra is defined by the following commutation g+ — ^/2( f^ f_, + f^ fr,) S^ ~ f^ f^ — f^ f-i 

relations , ^ , , , , 

P+ = ^/2(//^ - /t/_0, P' = -(/oV,s + flfo), 

[Sj,Sk] = iCjkiSi, [Pj,Pk] = iejkiSi, [Pj,Sk] = icjkiPi, ]\j+ = V2{f\fe - Pj-i), M' = -(/o/e + flfo)- 

A = i{flfs-flfe) (A9) 

[Mj,Mfe] = le.kiSi, [M,,Sk\ - le.kiMi, [P„ A4] = ^5^^, ^j^^^^ ^t^ ^t denote a creation operator of the fermion 

with spin "up" and "down" respectively whereas /o,/s, 
stands for spinless fermions. Then the Casimir operator 
[Pj , A] = iM, , [A, Mj] = tPj , [Sj ,A]=0. (A5) <EZl) transforms into the constraint 

(j, k, I are Cartesian coordinates). Besides, m + no + n-i + n, + Ur = 1 (AlO) 

S • P = 0, S • M = 0, (A6) 
I 



APPENDIX B: VERTEX CORRECTIONS AND WARD IDENTITIES 

We start with the classification of the vertex corrections containing one fluctuator line. The leading parquet 
diagrams are plotted on FigEla-i. These diagrams belong to three different topological classes drawn on lines 1-3 of 
Figl^l Following the standard Feynman codex we write the expressions for diagrams Figl^a-c 

TTT{u^,e)^{aj2fjfT^ J2 ^^^iS^' E E ^/^^^ (^-^ + ^i)G/3V^('^-^ + ^2-e3)g.,(ki, 6i).g.,(k2, e2)i?(e3) 

criCT2MlM2 eie2e3kik2 

(/3 ~ T,S). This expression corresponds to Fig|B^, other expressions for diagrams b-c belong to the same topolog- 
ical structure and differ by indices only. The summation is taken with respect to fermionic Matsubara frequencies 
£1,62 = 2z7rT(n -I- 1/2) and bosonic frequency €3 = 2iTrTn. Here the bare vertex J^^, which enters the diagrams 
(c,f,i), corresponds to the term ^'^'^X^fcfc' '^a-^^^'^ka-'^k'a ™ ^^"^ potential scattering term HpJ omitted in the SW 
Hamiltonian H24() . The tensor RZa[Zl'li' stands for kinematic factors containing S and P operators in a scalar product 
with Pauli matrices a and is given by, for example, 

^J^Misp'^ = i^^n^i ■ 0■crla){'Ps^n ■ (T„^^J(Pf^,s ■ cr^/crj, (Bl) 

-'^fj.fj.isjj.' ^ (S/JlM ' '''o-lo-2)(PsA'l ' ""o-'crJ^P/j's • CT 0-2(7) (B3) 
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T' T S ^T T/ S T \T T^ S S ^T 

(a) (b) (c) 






FIG. 6: First leading parquet corrections to TT vertex. The mapping of TT vertices on TS vertices is discussed in the 
text. Similarly to Fig. 5 the absence of direction in electron propagator lines assumes that both directions (clockwise and 
anti-clockwise) are possible. 



for Fig.6.a,d,g respectively. We assume also that the electron is taken on the mass shell while e is an energy transfer. 
The energy of the triplet state is also assumed to be small oj <C ^st- 
For diagrams Fig|51d-f one obtains 

TTT{L0.e)^{aj2fjfT^ J2 KZ^' E E ^/^A^i ('^-^+^i)G/3>^ (^-^+^1-^2-63)3., (ki, 6i).g,,(k2, e2)i?(e3) 

o"lO'2MlA'2 «ie2e3kik2 

and for Fig^lg-i one gets 

rTTiu;,e)^iaj2fj!^^'T^ Y, Kl\li' E E ^/3A-i(^-^+^i-^3)G^'p.(^-e+ei-e2),9.,(ki, 6i)5<,,(k2, 62)i?(e3) 



^1^2MlM2 



eie2e3 kikg 



Applying the procedure of analytical continuation explained in Appendix C and taking the limit e ^ one gets the 
following estimate for the real part of the diagrams Fig|SJa,b,d,c,g,h) 



Re[porTT(w,0)] ~ (aJ2)^.7iln 



D 



max[r, ti>]y \A 



In 



D 



^ST 



while for diagrams Fig|HJc,f,i) 



Re[poTTT{^,0)] ~ (aJ2)^jsln^ 



D 

A~~ 



ST 



These corrections are however parametrically smaller than the main Kondo diagram Fig. ^: 



Re[por^T(w,0)]~(ji)-^ln- 



3l„2 



D 

max[r, tj] 



(B4) 



(B5) 



(B6) 
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under realistic conditions a < 1 and Tk ^ ^st- 

The imaginary part of all diagrams a-i has a threshold form for uj ^ T 

ImrTT(tc',0) ~ {aj2fjfFpf,,{\uj\ ~ ^sT,T)e{]Lo\ - Ats) 



(B7) 



where .FJg^' (|w| — Ast, T ^ lu) ^ {\uj\ — Ast)" , and i^ = 1, 2 (see Refi^ for details of calculations). This estimation is 
done with accuracy 0(exp (—Ast/T)). In contrast, the imaginary part for Kondo vertices is given in the limit w 3> T 



Im[por^T(^,0)]^(ji)3signH 



(B8) 



The topological structure of singlct-triplct vertices is the same as on Fig|Sl The estimation for these diagrams gives 
the following expressions 



Re[porTs(w,0)]-02)'ln 



D 



[T,c 



maxJ , uj 



In 



D 



^ST 



and 



Next to leading (non-parquet) diagrams are shown on Figd 



(B9) 



(BIO) 



\ y 



\ v 



\ V 



\ / 







(a) (b) (c) (d) 

FIG. 7: Next to leading (non-parquet) corrections to the vertices. 
The analytical expression for diagram FigQa (and similarly for b) is given by 

rTT{Lu,e)^iaJ2fj^''T^ Y, RZ'"^'' Y. E ^^(^ + ^2)65(62)3.1 (ki, €1)5., (ka, ei + £3 - c^)^(e3) 

a-icr2fJ.ifJ,2 eie2e3kik2 

and for Fig[7|c (and similarly for d) 

TTs{io,e)^iaj2fJ^'T' Y Kl\7" Y. E ^^^^^^^ + ^2)65(62)3.1 (ki, ei)5.,(k2, ei + e2 + eg - a;)i?(e3) 

cri<T2fiifJ.2 C1C2C3 kik2 

The easiest way to calculate diagram Fig[7|a is to use the Ward Identity. The corresponding estimation is given in 
Section IV. As one expected, the imaginary parts of these vertex corrections have also threshold form {\uj\ > Ast)- 



<b} 



2„■3^ 



T/A 



ST, ^ST 



Ast <T <A 



ET 



Mp^t's¥\ = {^'A 



ujjAsT, T < Ast 



The damping H/t in the frequencies interval \lo\ < Ast is exponentially small (~ cxp(— As^/T)) being proportional 
to the population of the singlet state. 

The real parts of TS vertices Fig[7|c-d are estimated as 



Re[porsT(cc',0)]-(a2j|)ln 



D 



iST 



(Bll) 



The perturbative results including analysis of leading log and sub-leading log diagrams summarized in this section 
allows to incorporate noise corrections associated with the vertex fluctuator into standard one- and two- loops RG 
approach. 
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APPENDIX C: ANALYTICAL CONTINUATION 

We start with the derivation of the general equation for the vertex corrections Tap{w, e) : 

r.,(c., e) = r'^'r'P'jP'P J g [lmG^,(z)G« (^ + e)n^(c. - z) tanh [^) 
G^, (z)ImGf, {z + e)n«(Lj - z) tanh {^^] + Gf, (z)G^, {z + e)lmR^{uj - z) coth 



2T 
Here the arguments in the vertex are complex variables defined in the upper half-plane, 

n(z) ^ \y ^o*^ [if] [^"^(^ " e)Im-D^(e) + D^\z - e)Imif^(e) 



2T 



(CI) 



(C2) 



Then the elastic diagonal vertex Vtt reads 



r^T(a;,0) = (j^^)^J^^|g 



(Gf(z)) Imn^(a;-z)coth 



2^ .2ReGf(^)ImGf(^)n«(a;-.itanh(^) 



This result may also be obtained by means of the Ward identity (|52(l for the self energy insert S5t(w) in the Green 
function Gt^ihi2-^ which enters the triplet vertex Fig. [3 



grj,^2 / de \-,^^rRf.\TiRf^ .\^^„i. ( '^ \ i T.^(~ 



^STiz) = {J^'Y I - ^ImL«(e)i?^^(z-6)tanh^2ry 



L«(z-6)ImI?«(e)coth(^) 



where 



L^iz) 



de 

2^ 



ImG«(e)7^^(z - e) tanh ( ^) + G^{z ~ e)lmK^{e) coth (^) 



V2Ty 



(C3) 



(C4) 



The imaginary part of Eq. (|Cip gives the expression (|48|l for the dephasing rate 
No Ward identity exists for the non-diagonal vertex 

Trsico, e) = {j'^f J g [imGf (z)G^(. + e)n^{co ~ z) tanh (^) 
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